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We show that a difference in time scales of vertical and lateral dynamics permits one to analyze 
the problem of interacting electrons confined in an axially symmetric three-dimensional potential 
with a lateral oscillator confinement by means of the effective two-dimensional Hamiltonian with a 
screened Coulomb interaction. Using an adiabatic approximation based on action-angle variables, 
we present solutions for the effective charge of the Coulomb interaction (screening) for a vertical 
confinement potential simulated by parabolic, square, and triangular wells. While for the parabolic 
potential the solution for the effective charge is given in a closed anlytical form, for the other cases 
similar solutions can be easily calculated numerically. 
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I. INTRODUCTION 

It is well known that there is a restricted class of ex- 
actly solvable problems in quantum mechanics. Such 
examples serve as paradigms to illustrate fundamental 
principles or/and new methods in their respective fields. 
This is especially important for finite systems, where ap- 
proximative methods are indispensable to a many-body 
problem. In particular, two-electron systems play an im- 
portant role in understanding of electron correlation ef- 
fects because their eigenstates can be obtained very ac- 
curately, or in some cases, exactly. 

The most popular model to study the electronic 
exchange-correlation energy in a density functional the- 
ory is the Hookean two-electron atom (HA). The basic 
HA is two electrons interacting by the Coulomb poten- 
tial but bound to a nucleus by a harmonic potential that 
mimics a nuclear-electron attraction. For certain values 
of the confinement strength, there exist exact solutions 
for the HA ground state [l|, Q . When the HA is placed 
in a perpendicular magnetic field. Taut provided analyt- 
ical solutions for a two-dimensional (2D) HA at partic- 
ular values of the magnetic field [3[. This model can be 
equally viewed as a model 2D quantum dot (QD). Recent 
progress in semiconductor technology made it possible 
to fabricate and probe such confined system at different 
values of the magnetic field 0, Q • Consequently, it has 
stimulated numerous theoretical studies on two-electron 
QDs, so-called "artificial atoms" (see for a recent review 
16]). Being a simple nontrivial system, QD He poses a sig- 
nificant challenge to theorists. For example, using a 2D 
He QD model, one is able to reproduce a general trend 
for the first singlet-triplet (ST) transitions observed in 
two-electron QDs under a perpendicular magnetic field. 
However, the experimental positions of the ST transition 
points are systematically higher 0, 0]. The ignorance 
of the third dimension is the most evident source of the 
disagreement, especially, in vertical QDs [1, H, [T^l- 



Although accurate numerical results for QD He can 
be obtained readily, analytical results are still sought, 
because they provide the physical insight into numeri- 
cal calculations. Moreover, analytical results could es- 
tablish a theoretical framework for accurate analysis of 
confined many-electron systems, where the exact treat- 
ment of the three-dimensional (3D) case becomes compu- 
tationally intractable. The purpose of the present paper 
is to introduce a consistent approach which enables one 
to reduce the 3D Coulomb problem to the effective 2D 
one, without loosing major effects related to the QD's 
thickness. In Section II we develop the concept of the 
effective charge of the Coulomb potential based on the 
adiabatic approximation [ll| for axially symmetric 3D 
systems. The calculations of the effective charge for a 
Hamiltonian with a vertical confinement approximated 
by parabolic, square, and triangular well potentials are 
presented in Section HI. Finally, we will show in Section 
IV that recent experimental data ^ can be successfully 
reproduced within a 2D approximation but with addi- 
tionally screened Coulomb interaction due to the thick- 
ness of the sample. The main outcomes are summarized 
in Section V. Three Appendixes provide some technical 
details of the calculations. 



II. BASIC REMARKS 
A. Model 

The system Hamiltonian for the case of a magnetic 
field B along a symmetry axis z reads 
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Here the term Vc = k/\ri—r2\ with k = e'^/Aireoer 
describes the Coulomb repulsion between electrons and 
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Hspin = 5*Mb(si + S2)B is the Zeeman term, where 
jJiB = \e\Ti/2m(.c is the Bohr magneton. For the mag- 
netic field we choose the vector potential with gauge 
Ai — iB X Yi — ^B{—yi,Xi,0). The confining poten- 
tial is approximated by a 2D circular harmonic oscilla- 
tor (HO) in xy-plane and the vertical confinement Vz'- 
U{v,) = m*ujl pi/2 + Vziz,); rf = pi + zf, pi = xi + yf 
and hujQ is the energy scale of confinement in the xy- 
plane. Below we analyze different forms for the vertical 
confinement Vz- 

For our analysis it is convenient to use cylindrical coor- 
dinates {p,4>,z). Also, we separate the Hamiltonian ^ 
on several terms: H = Hq + Hz + Vc + Hspin, where 
the term Hq = consists of the contributions re- 

lated only to the lateral dynamics (x?/— plane) of non- 
interacting electrons 
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Here, the effective lateral confinement frequency il 
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^)^/^ depends on the magnetic field by means of 



the Larmor frequency ujl = \e\B /2m*c\ Iz^ = Ptpi is the z- 
component of the angular momentum of the i-th electron. 
The eigenstates of the single-particle Hamiltonian ((2]) arc 
well-known Fock-Darwin states (cf [12j). The motion of 
noninteracting electrons in the z-direction is described by 
the Hamiltonian Hz = Z]i(P.,/2m* + Vz{zi)). Since the 
magnetic field is directed along the z axis, the Zeeman 
term is Hspin = 9* P'bSzB . This term is not important 
for our analytical study and will be taken into account 
only in numerical analysis of experimental data. 



B. Adiabatic approximation 

For typical QDs {Tiloq ~ 3 meV) the contribution of the 
Coulomb interaction to the total energy is comparable to 
the confinement energy at zero magnetic field Q. Evi- 
dently, the standard perturbation theory is not valid in 
this case. In real samples the confining potential in the z- 
direction is much stronger than in the ccy-plane. It results 
in different time scales (see below) and this allows one to 
use the adiabatic approach [ll[ . To lowest order the adi- 
abatic approach consists of averaging the full 3D Hamil- 
tonian over the angle variables 6*.. = w^.t (fast variables) 
of the unperturbed motion (fc — 0) after rewriting the 
{zi,Pzi) variables in terms of the action- angle variables 
{Jzi,Ozi)- As a result, the motion effectively decouples 
into an unperturbed motion in the vertical direction gov- 
erned by the potential J2i ^{Jzi, OzJ and into the lateral 
motion governed by the effective potential 

VM{x, y}; {Jz}) ^^v, + V.^{p- Jz,,Jz.), (3) 



where v is defined in Eq.([2]), p = [{xi 
y2?Y'^ and 
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is the effective electron-electron interaction that contains 
the memory on z dynamics through integrals of motion 
Jz- ■ The effective interaction affects, therefore, only the 
dynamics in the lateral plane, where the confining poten- 
tial is the parabolic one (see Eq.Q). Hence, the effective 
Hamiltonian for two-electron QD reads as 



HcS — Hq + Ez + V:° 



(6) 



(5) 

where Ez — '^^St and Si is the electron energy of the 
unperturbed motion in the vertical direction. 

Our ansatz consists in the consideration of the effective 
interaction (jH) in the form Vj^f = kf{p)/p. Then one can 
define the effective 2D Coulomb interaction 

fccff 

P 

where the effective charge is the mean value of the factor 
f{p) upon the nonperturbed lateral wave functions, i.e., 

hs = k{f{p))^{pV^«ip)). (7) 

In contrast to a standard 2D consideration of the bare 
Coulomb potential {Vc — k/p) in QDs d, @], in our ap- 
proach the electron dynamics is governed by the same po- 
tential but with the additional effective (screened) charge 
due to the QD's thickness. One of the main advantages 
of this approach is that the interaction matrix elements 
can be expressed in an analytical form. Thus, we shall 
use (diagonalize) the effective Hamiltonian ^ with the 
effective Coulomb interaction V^^ , i.e., =J> V^^ . 

According to the Kohn theorem the center of mass 
(CM) and the relative motion of the 2D system described 
by the Hamiltonian ([5]) are separated and the mean value 
in Eq.© can be evaluated using the Fock-Darwin states 
for the relative motion 
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This state is eigenfunction of the operator with 
eigenvalue m and the radius-dependent function with a 
radial quantum number Up has the form 
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where p = m*/2 is the reduced mass, f — {fj,fl/h)-'-/^p 

and il™' denotes the Laguerre polynomials [l^. For the 
lowest states (with different values of the quantum num- 
ber m but with the radial quantum number Hp = 0) one 
obtains for the effective charge 



2 (pVL 



|ni|+l noa 



(10) 



3 




FIG. 1: Left: the localization of QD in the layer of the thick- 
ness a. Right: the schematic representation of the position 
of zero-point motion in the parabolic confinement relative to 
the layer thickness. 



III. EFFECTIVE CHARGE 

A. Parabolic potential 

To account for effect of localization of the dot in the 
layer of thickness a, let us first consider a 3D model with 
a vertical confinement approximated by a parabolic po- 
tential Vzizj) — i m*LL!l zf (see Fig.[T]). Due to the Kohn 
theorem [1J| the CM and the relative motions are sep- 
arated. The solution for the CM motion is well known 
(cf [13 )■ It is not important for our discussion, since 
the CM dynamics does not affect the electron interac- 
tion. The 3D Hamiltonian for the relative motion of two 
electrons has the form 
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The term /irci is defined by Eq.(l2|) in which the effec- 
tive electron mass is replaced by the reduced mass /x; all 
single-electron variables (with index i) are replaced by 
the corresponding variables for relative motion (without 
indices). 

After rewriting the {z,pz) variables in terms of the 
action-angle variables {Jz,dz) 
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we integrate out of the fast variable, i.e., average the 
Hamiltonian (fTT|) over the angle dz- As a result, the 
effective interaction potential (sec also Appendix [A| is 
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where K{x) = /J^^^(l — x sm"^ 6)~^^'^d9 is the complete 
elliptic integral of the first kind (see flBl)- This integral 
is well defined for all values of p. 

The effective Hamiltonian for the relative motion is 
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To replace the effective electron-electron potential V^'^ 
by the Coulomb- type Vq^ in Eq. fM)) . we must determine 
the effective charge. Taking into account the definitions 
Eqs.lHl), (113 and the result Eq.lUS]), one obtains the fol- 
lowing for the effective charge 
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where \np,m) are the Fock-Darwin states for the relative 
motion, Eqs. dS]),®. 

For the lowest states (up = Uz — 0) one reduces 
Eq. p^ to the form of integral (fTU]) . As a result, the 
effective charge can be expressed in terms of the Meijer 
G-function ,16.] 
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Guided by the adiabatic approach, it is instructive 
to compute the effective charge by dint of quantum- 
mechanical mean value of the Coulomb term in the 3D 
oscillator state \np,m)\nz) 

= {{pVcip, z))) = k (((1 + zyp')-^/')) . (18) 

Here, \nz) is a normalized one-dimensional harmonic os- 
cillator wave function p^ . Since the lateral extension 
exceeds the thickness of the QDs by several times, one 
may suggest to consider the ratio {z/ pY as a small pa- 
rameter of theory. Note, however, that the averaging over 
the 3D oscillator state jrip, to) \nz) implies the application 
of the first order perturbation theory for calculation of 
the contribution of the Coulomb interaction in QDs. For 
Ur, — n, — Q one obtains 
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where Kq is the modified Bessel function of the 2nd fective charge Eas. (fT6)l . p9|) there is a contribution of the 
kind. One observes that in both definitions of the ef- electron dynamics along the coordinate z. Below we will 
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compare true contagion of both definitions upon the in- field < fJ/Wz < 1, available in experiment, the solution 
terpretation of the experimental data. for the last integral can be expressed in terms of the 

For small and relatively large values of the magnetic hypergeometric functions 
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B. Infinite square well 



C. Triangular well 



A hard wall potential is among popular models for the 
confinement in QDs. In case of the vertical confinement, 
Vz is the one-dimensional infinite square well 



0, \z\<a/2 
oo, \z\ > a/2 



(21) 



where a is the thickness of the layer in which the dot is 
created (see Fig. [T]). 

Due to uniform motion inside the square well, the an- 
gle variable 0^. is proportional to the coordinates Zi (see 
Eq. (jBSp in Appendix [B|) and, therefore. 
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Integrating over the coordinates zi and Z2 between the 
walls of the potential (I^T]) one obtains the effective 
Coulomb interaction V^*^ — kA{£,)/a, where ^ = p/a 
and 
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The effective Hamiltonian now reads 
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where the term hi is determined by Eq.Q and we use 
Eq. (|B2p for the contribution of the vertical confinement. 
In contrast to the parabolic potential, here does not 
depend on the quantum numbers riz- . In order to re- 
place — > V^^ we must define the effective charge. 
Similar to the previous case, the effective charge is the 
mean value (pl^J;f(p)) in the Fock-Darwin states. For 
the lowest states with different m we obtain 
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where b = pfla^ /h. 
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Let us consider the well potential characterized by an 
infinitely high barrier for z < and a linear potential 
Vz{z) = eFz for z > 0; the product of the electron charge 
e and an electric field F is assumed to be positive: 
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This potential is a simple realistic description of the po- 
tential well at doped heterojunction or/and in the case 
when an external voltage is applied to the (top and bot- 
tom) electrodes of the sample in which the dot is created. 

By means of the relation Eq. (|C6|) (see Appendix [C| 
between the coordinates Zi and the angle variables Oz^ , 
for the case when both electrons occupy the same energy 
level e, the effective Coulomb term is 
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where ^1 = 01-6*2, C2 = 6*1+02-2, 6 = p/zm, z„r 
and 9i = OzJtt. 

According to Ea. (|C7p . the lowest {uz = 1) energy level 
in the triangular potential is ei = ci{heF/^/2m*y^^ , 
where ci = (Qtt/S)^/'^ ~ 2.32. The parameters can be 
selected to ensure the corresponding wave function (the 
Airy function, the first zero of this function occurs at 
ci — 2.338; see, for example, ReflTsI) will be extended 
until the desired length in the vertical direction. Thus, 
in the adiabatic approximation the effective Hamiltonian 
has the form 



Ho + 2ei + V^,^{p) 



(28) 



where the term Hq is defined in subsection III Al As 
above, we replace 
charge. For the states with Up 



Vq^ and calculate the effective 
(the Fock-Darwin 
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states for the relative motion) we obtain the effective 
charge (see Eq.([71)) 

fceff = fci^ / d6 e-^«3 4^l™l+^6(6) , (29) 
Jo 

where b = fAlzm^ /T% and Zm = £i/eF. We recall that the 
parameter is the distance between turning points. It 
is similar to the parameter a in the square well potential. 
The integral over the variable ^2 in the expression for B 
fEg. ipT)) ) can be evaluated analytically but the other two 
in Eqs. (P7|) and must be calculated numerically. 



IV. DISCUSSION OF RESULTS 

In Fig. [5] we present results for the effective charge (ra- 
tio keff/k), calculated for the lowest states (rip = 712 = 0) 
with different to, with the aid of: alEq. pT)) for the 
parabolic confinement as a function of the ratio Q/cOz 
(see Fig.[i;a)); b)Eq.(Iin]) and Eq.d^n]) for the square and 
triangular well potentials, respectively, as a function of 
the variable b (see Fig. [H^b)). These results evidently 
demonstrate that the inclusion of the vertical dynamics 
reduces the Coulomb interaction. In all considered mod- 
els for the vertical confinement this effect affects strongly 
quantum states with small values of the quantum number 
TO. Such states are major participants in the ground state 
transitions at small and intermediate values of the mag- 
netic field. Therefore, the attenuation of the Coulomb 
interaction due to the sample thickness may explain the 
disagreement between experimental data and predictions 
upon the position of the singlet-triplet transitions based 
on calculations with the 2D Coulomb potential (see be- 
low). In other words, the 2D calculations (with the 
charge k) overestimate the electron correlations in two- 
electron QDs. 

Although the results for the effective charge in the 
parabolic potential (adiabatic approach) (see Fig[2^) are 
slightly different in comparison with those of the hard 
potentials, a general trend for this charge determined by 
the magnetic field is similar (compare Figsl2^,b). Notice 
that the magnitude of the ratio and its evolution with 
the increase of the magnetic field are almost the same 
in the square and triangular potentials (see Fig[2j3). In 
contrast to the plain quantum-mechanical averaging pro- 
cedure (sec Eqs.([T^- l^ ) in the adiabatic approach (see 
Eq. lfTT)) ) the Coulomb interaction (the ratio kes/k) is re- 
duced stronger (see below). 

To elucidate the quality of the adiabatic approach we 
use the 3D axially symmetric oscillator model with the 
lateral {hujQ — 3meV) and the vertical {huj^ — 12meV) 
confinements and compare results for energy levels in 2D, 
3D and the effective charge (adiabatic) approximations. 
For a numerical analysis of spectral properties we choose 
the effective mass to* = 0.067TOe, the relative dielectric 
constant of a semiconductor Sr = 12 and \g*\ — 0.44 
(bulk GaAs values). In all three approximations we solve 



the eigenvalue problem by means of exact diagonalization 
of the corresponding Hamiltonian in the Fock-Darwin ba- 
sis. We recall that in this basis the Coulomb matrix 
elements are given in analytical form (cf _12,j). One ob- 
serves a remarkable agreement between results for the 3D 
and the effective charge approximations (see Fig. [31). In 
the 2D approximation the Coulomb interaction is much 
stronger and, consequently, the evolution of levels with 
the magnetic field as well as their absolute values are 
different from those of the 3D approximation. 

To illuminate the key advantage of the effective charge 
concept it is noteworthy to analyze the available exper- 
imental data Q within various approaches. To this aim 
we will compare the results of calculations of the addi- 
tional energy A/i = fi{N) - n{N - I) [3, l5| in the 2D 
approximation (with the effective Coulomb interaction) 
and in the 3D approach with a full Coulomb interac- 
tion. Here, the chemical potential of the dot /i(iV) = 
E{N) — E{N — I) is given by the ground-state energy of 
the dot E{N) with N and iV — I electrons. To conduct 
this study, one needs to estimate the vertical confinement 
frequency fiWz in a real sample. The simplest approach 
is to use the parabolic model for the vertical confinement 
in the layer of thickness a (see Fig[T]) . 

The thickness of QDs is much smaller in comparison 
with the lateral extension. Therefore, the vertical con- 
finement hujz is much stronger than the lateral confine- 
ment TiLj{) and this fact is, usually, used to justify a 2D 
approach for the study of QDs. However, there is a 
nonzero contribution from the vertical dynamics, since 
the energy level available for each of the noninteracting 
electrons in the z direction is e = hu)z{nz + 1/2). For the 
lowest state nz — ^ ei ~ ^ hujz- By dint of the condi- 
tion Vz{±Zm) = '7i*'^z^m/2 = El onc dcfincs the turning 
points: Zm — (1^1*^1 z) (see FigH]). We assume that 
the distance between turning points should not exceed 
the layer thickness, i.e., 2zm < a. From this inequality it 
follows that the lowest limit for the vertical confinement 
in the layer of thickness a is 



or b > 2il/u!z- For typical GaAs samples with the thick- 
ness a between 10 nm and 20 nm this estimation gives the 
minimal value for huJz between 45meV and If meV, re- 
spectively. These estimations provide a genuine cause for 
the use of the adiabatic approach in case of QDs, since 
T,(= 2tt/lUz)<^To{= 2n/LJo). 

Using the " experimental" values for the lateral confine- 
ment and the confinement frequency ujz as a free param- 
eter, we reproduced successfully with the value hojz — 8 
meV and \g*\ = 0.3 the positions of kinks in the addi- 
tional energy for a two-electron QD (see Fig|4^) 

Afi = E,,i - E{1) + E.pin (31) 

in all three sarnples in the 3D axially symmetric os- 
cillator model [10l|. Here, E'roi = (-ffrci) is the relative 
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FIG. 2: (Color online) The ratio k^s/k for the lowest states with different m for the vertical confinement approximated by: 
(a) the parabolic potential; (b) the hard wall (square/triangular well) potentials, (a) The results obtained in the adiabatic 
approximation and by means of the plain quantum-mechanical averaging procedure (see subsection IIII A[) are connected by 
dashed (blue) and dotted-dashed (green) lines, respectively, (b) The results for the square and triangular well are connected 
by solid (pink) and dotted-dashed lines, respectively. 
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FIG. 3: Lowest energy levels (only the relative motion is con- 
sidered) of the two-electron QD with fully parabolic potential. 
The results of the 2D, the 3D and the effective charge (adia- 
batic) approximations are connected by solid (orange), solid 
(red), and dashed (blue) lines, respectively. The result of 
the plain quantum-mechanical averaging for the m = level 
(fceff = {{pVc))) is connected by the dotted-dashed (green) 
line. 



energy, ^^(l) = huJo + hujz/2 is a single-electron energy 
and the Zeeman energy Espin = -\g*\nB[^ - {-l)™']B/2 
is zero for the singlet states. Note that it was found 
from the Zeeman splitting at high magnetic field that 
|Sf*| = 0.3 While at small magnetic field the re- 

sults for \g*\ = 0.44 and \g* \ = 0.3 are similar, the latter 
value provides the best agreement between the full 3D 



calculations and the experimental data. 

We recall that in the 2D approach used by Nishi et al. 
one encounters the problem of the correct interpre- 
tation of the experimental data (see FigH^, a plain lat- 
eral confinement hujQ = 2.9 meV). In the effective charge 
approximation the vertical confinement is taken into ac- 
count with the aid of keS in the 2D effective Hamiltonian. 
The remarkable accord between the predictions based on 
the results of Section IIII Al and the observation confirms 
the vahdity of the suggested concept (see Fig. |4]). We 
stress that the effective charge approximation facilitates 
the calculations providing the same results as the full 3D 
calculations. 

Notice that the results based upon the adiabatic ap- 
proximation are in a better agreement with the full 3D 
calculations in contrast to those obtained with the aid 
of the plain quantum-mechanical averaging procedure 
for fccff = {{pVc{p, z))). As discussed in Section IlIBI 
the adiabatic approach is based on the effective separa- 
tion of fast (vertical) and slow (lateral) dynamics with 
subsequent averaging procedure. In contrast, the plain 
quantum-mechanical averaging represents a type of per- 
turbation theory based upon the first order contribution 
with respect to the ratio z/p only. The higher order term 
may improve the agreement at small magnetic field, since 
the vertical dynamics is non negligible and affects the 
lateral dynamics. The increase of the quantum number 
TO, caused by the increase of the magnetic field strength, 
reduces the orbital motion of electrons in the vertical di- 
rection. The larger is m the stronger is the centrifugal 
force, which induces the electron localization in a plane, 
and, therefore, the lesser is importance of the vertical 
electron dynamics. In the limit of strong magnetic field 
(large m) the dot becomes more of a "two-dimensional" 
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FIG. 4; (Color online) (a) The additional energy Afi as a func- 
tion of the magnetic field in the parabolic model for a two- 
electron QD. The results of calculations with a lateral con- 
finement only (the 2D approach, Tiuq = 2.9 meV, \g*\ = 0.3) 
and full 3D approach [10| {huJz = 8meV) are connected by 
solid (orange) and (red) lines, respectively. The vertical grey 
lines indicate the position of the experimental crossings be- 
tween different ground states f7l|. Ground states are labeled 
by {m,S), where m and S are the quantum numbers of the 
operators Iz and the total spin, respectively. The results 
based upon the adiabatic approximation, Ea. (|17|l . and the 
plain quantum-mechanical averaging procedure, Ea. (|19|l . are 
connected by dashed (blue) and dot-dashed (green) lines, re- 
spectively, (b) The ratio kas/k as functions of the magnetic 
field based on Eq. (|17p and the plain quantum-mechanical av- 
eraging, Eq. (fT9)) are connected by dashed (blue) and dotted- 
dashed (green) lines, respectively. 



system. This explains the improvement of the accuracy 
of the plain quantum-mechanical averaging procedure at 
large m, i.e., for the ground states at high magnetic fields. 

The results of different averaging procedures indicate 
that the accuracy of the methods depends, indeed, on 
the quality of the separation between the lateral (p) and 
vertical (z) components of the Coulomb interaction, with 
further averaging over them. We recall that the essence 
of the "removal of resonances" technique based on the 
action-angle variables [ll| is the identification of adia- 
batic invariants, i.e., approximately conserved integrals 
of motion. The averaging of the Hamiltonian over the 
fast variables is a next step within this technique. In our 
case such an invariant is the action J^, which is approx- 
imately conserved due to large difference in a time scale 
between fast vertical (z) and slow lateral (xy) dynam- 
ics. Therefore, the integration over fast variable (angle) 
is well justified in this case. As it was pointed out above, 
at small magnetic field a plain quantum-mechanical pro- 
cedure, which reduces to a variant of the first order per- 
turbation theory, is less justified. In fact, the effective 
potential obtained from the plain quantum-mechanical 



FIG. 5: (Color online) (a) The magnetic dependence of the 
additional energy Afi for a two-electron QD with the lateral 
confinement Hljo = 2.9 meV and \g*\ = 0.3. The results of 
the pure 2D approach are connected by the solid (orange) 
line. The results of the adiabatic approach (effective charge 
calculations) are connected: by the dashed (blue) line for the 
parabolic model {huiz = 8meV); by the solid (pink) line for 
the square well model (a — 32.5 nm); by the dotted-dashed 
(black) line for the triangular well model {F — 220kV/m). 
See the text for details. The vertical grey lines indicate the po- 
sition of the experimental crossings between different ground 
states labeled by (m, S) (b) The ratio kes/k for three dif- 
ferent approximations of the vertical confinement as a func- 
tion of the magnetic field. 



procedure is significantly stronger at small distances be- 
tween two electrons (see Appendix E)) . In turn, in this 
case the effective coupling between the vertical and lat- 
eral confinements is stronger, which indicates a less effi- 
ciency of the separation procedure. 

The positions of the singlet-triplet transitions in the 
same sample can be reproduced as well with a vertical 
confinement approximated by the square (Section IIIIBp 
and the triangular (Section IIII Cp potentials. It appears 
from the behavior of the effective charge for the consid- 
ered models of the vertical confinement (see Fig. ^ that 
one obtains equivalent results if the ratio k = b/{n/ujz) 
between the distance b in the square/triangular well and 
the parabolic variable ^l/ujz is k ~ 4. Note that this 
conclusion is independent on the particular value of the 
vertical confinement in the parabolic model, since this 
parameter enters in the calculations of the fccfr through 
the ratio Q/ujz- 

Within the square well model the best agreement with 
the experimental data is obtained for a « 32.5 nm (see 
FiglS]). In virtue of the definition b ~ ^i^la? /h we have 



n = bl[- 



2r 



Using the above value for a and TilJz 



(32) 



! meV which 
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provides the equivalent result in the parabolic model, we 
obtain k = 3.715. Indeed, this estimation is close to the 
one extracted from FigO which is independent on the 
specific parameters of the QD. 

Similar agreement with the experimental data is ob- 
tained by means of the triangular well potential (6 = 
fi^lz'^/h) for the electric field F « 220 kV/m which corre- 
sponds to Zm ~ 32.1 nm (see Fig. [5]). Note that in all con- 
sidered models the effective charge k^s of the Coulomb in- 
teraction changes with the increase of the magnetic field 
in a similar way (see Fig. [51) producing approximately the 
equivalent effect in all three cases. At very large magnetic 
fields the effective charge approaches almost the value 
fceff « 0.95fc. 

Finally, a remark is in order. For the sake of il- 
lustration we have used only fcofi for the lowest basis 
states (rip — n'^ = 0) with different m. However, even 
for the ground state calculations with the aid of the 
exact diagonalization, the interaction matrix elements 
kc{f{n'p,m\p^^\n'p,m) with np^n'^ > have been taken 
into account. Namely, the diagonalization is performed 
using the interaction matrix elements up to np,n'p = 10. 
Obviously, the accuracy of the method would improve 
if one calculates the 'effective charge matrix elements' 
~ (^pj ™IP^nt'l"'pi ™) for each interaction matrix 
element. However, in this case the procedure would lose 
the simplicity and become impractical. Fortunately, from 
the comparison of the present results with exact 3D cal- 
culations we found that for the analysis of the ground 
state properties it is sufficient to use only the elements 
^Q™'' = even for the interaction matrix elements 

with Up, n'p > 0. 



V. SUMMARY 

We developed the effective charge approach taking full 
account of the thickness of two-electron quantum dots. 
Our approach is based on the adiabatic approximation 
where the full 3D dynamics of two interacting electrons 
is separated by means of action-angle variables on the in- 
dependent vertical motion and the lateral dynamics de- 
scribed by the effective 2D Hamiltonian. The separation 
is reached due to different time scales in the vertical (fast) 
and lateral (slow) dynamics and it is well justified in all 
types of QDs (vertical and lateral). As a result, one must 
solve only the Schrodinger equation for the 2D effective 
Hamiltonian where the full charge k is replaced by fccg 
(see Eqs.©, (flUl) '). The eigenvalue problem was solved 
by means of the exact diagonalization of the effective 2D 
Hamiltonian in the Fock-Darwin basis. To demonstrate 
the feasibility of the effective charge approach, we consid- 
ered three different models for the vertical confinement: 
the parabolic, the square and triangular well potentials. 
The use of the parabolic potential simplifies the analysis 
due to the Kohn theorem which results in the analyti- 
cal expression for the effective charge, Ea. pT)) . For the 



square and triangular well potentials we obtained expres- 
sions for the effective charge Eas. (P5l) . (P^ . respectively, 
that can be easily calculated numerically. 

The value of the effective charge depends on the (good) 
quantum number m of the correlated state (see Figj^J. 
We established a scaling factor, Eg. ([5^ . between the 
variables used in the parabolic and the square/triangular 
well potentials. Taking into account this factor, we found 
that the effective charge (screening of the Coulomb in- 
teraction) affects quantum states for all potentials in a 
similar way. The screening due to the sample thickness 
is especially strong for quantum states with small val- 
ues of the quantum number m. We recall that these 
states determine the structure of the ground state tran- 
sitions at small and intermediate values of the magnetic 
field. Therefore, the screening provides a consistent way 
to deal with the effect of the thickness upon the posi- 
tion of the singlet-triplet transitions. In particular, the 
screening should be taken into account for the analysis 
of evolution of the energy difference between singlet and 
triplet states in the magnetic field. This energy is con- 
sidered to be important for analysis of the entanglement 
and concurrence in QDs (cf pigj). 

The comparison of the results with available exper- 
imental data 01 demonstrates a remarkable agreement 
and lends support to the validity of the approach. Be- 
ing important for the states with small quantum number 
TO, the screening of the Coulomb interaction becomes 
small for the states with large m, which dominate in 
the low-lying spectrum at large magnetic fields. On the 
other hand, these states cause strong centrifugal forces 
which induce the electron localization in a plane. In 
turn, the stronger is the magnetic field the less important 
is the vertical confinement. It follows that the 2D bare 
Coulomb potential becomes reliable in 2D approaches for 
the analysis of the ground state evolution of QDs only at 
very large magnetic fields. 

Thus, the 2D calculations with a bare Coulomb inter- 
action (with the charge k) overestimate the electron cor- 
relations in two-electron QDs at small and intermediate 
values of the magnetic field. One may induce that sim- 
ilar conclusion could be valid for QDs with more than 
two electrons. Notice, however, that the confining fre- 
quency in the lateral plane decreases with the increase of 
the electron number > 2 in exact 3D calculations at 
fixed vertical confinement [l^], making QDs to be more 
of a two-dimensional system. A thorough analysis of our 
approach for A^ > 2 is left for the future. 
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APPENDIX A: EFFECTIVE INTERACTION 
POTENTIAL FOR THE PARABOLIC VERTICAL 
CONFINEMENT AT = 

The effective electron-electron interaction potential is 
obtained by averaging the Coulomb term Vc{p, z) in the 
z-direction 



V^^{p)= I w{z)Vc{p,z)dz, 



(Al) 



where w{z) is the corresponding weight function. By 
dint of the quantum-mechanical averaging one obtains 
the effective potential 



cff(qm)/ N _ 



(P)= / \i'nA^rVcip,z)dz, 



(A2) 



while in the adiabatic approach we have the effective po- 
tential 



eff(ad) 



(P) 



1 



ir/2 



Vc{p,z{0,))d0, 



For riz — 0, we have for the quantum- mechanical weight 
function 



V^zo 



(A4) 



while, by virtue of Eq. (fT2|) and Jz — h/2, one obtains for 
the adiabatic weight function 



Wad = TT ^{d9z/dz) 



TTZoy/l- (z/zoY 



(A5) 



where zq = sjhj [luj^ ■ As it can be seen from Fig. [D^a) the 
weight functions (solid line) and Wad (dashed line) 
are completely different and, therefore, one may expect 
different contributions of the z-motion in the effective 2D 
potential, i.e., a different p-z coupling. 

By dint of the definition Vc = fc/(p^ + 2^)1/2 and 
Eas. (jA4p , (jA5P one obtains for the effective interaction 
potentials (Eqs. (jA2l) .((A3 |) ): 



ff(qm) 
it 

cff(ad) 



/TTZo 



ePV2.o^^^(ii_)^ (A6) 



Izt, 



■np \ p-" 



(A7) 



where Kq and K are the the modified Bessel function of 
the 2nd kind and the complete elliptic integral of the first 
kind, respectively. 

Although the potential V'^*^''^™-' is akin to the poten- 

tial the ratio R = V^i^f "'^Vv^f < 1 at 

< p < \/2zo (see Fig. Hfb)). We recall that for the 
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FIG. 6: (Color online) (a) The weight functions Wqm (solid 

line). Wad (dashed line) and (b) the ratio R — VJJj^'^'^V^inf'''™'' 
for the parabolic vertical confinement. 



lowest limit of the sample thickness a — 2zm — V^zq 
(see Figll]). In other words, when electrons are close to 
each other (for p < a), the effective quantum-mechanical 
potential is significantly stronger than the adiabatic one. 
At distances between electrons p > a two potentials reach 
each other asymptotically {R ^ 1). 



APPENDIX B: ACTION-ANGLE VARIABLES IN 
THE ONE-DIMENSIONAL SQUARE WELL 
POTENTIAL 

If £ is the energy of a particle of the mass m* mov- 
ing in the infinite square well potential (|2ip . its momen- 
tum pz = ±^/2m*e change the sign at the turning points 
±a/2. In this case the action variable is 



1 /""/^ a r- 

Jz — — Pzdz — —\f2m'e . 



a/2 



(Bl) 



Accordingly, the energy expressed in terms of is e = 
tt"^ J1 / {2m* a^) . By dint of the semiclassical quantization 
condition Jz = h{nz + 1) (n^ = 0,1,2,...), one obtains 
energy levels for a particle moving in the potential (I2ip . 
The value 1 in this condition is a Maslov index (cf [17|). 
It is convenient to define Hz ^ Uz -\- I. Finally, one has 



2m* 



TTTlz 

a 



uz = 1,2,, 



(B2) 
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Integrating equation m* z = pz, from the condition 
z{T/A) = a/2 (for -a/2 < z < a/2) one can define the 
period T = 2ayjm* /2£. It follows that the angle variable 
0, has the form 



27r vr / 2£ tt . , 

^ujzt^—t^-J—t^-z. (B3) 
1 a V in* a 



and, consequently, the period T = 2-k jiOz = 2\/ 2m* e / eF . 
Integrating equation m* z ~ pz, one obtains 



^ , dz \j2vn*e / eF , 



for < 2 < 



APPENDIX C: ACTION-ANGLE VARIABLES IN 
A TRIANGULAR WELL POTENTIAL 

For a particle moving with energy e in the triangular 
well potential ([26|) the turning points are z = and z = 
Zm- In virtue of equation pz = ±^j2m*{e — eFz) = 
one defines the variable Zm — s/ eF and, respectively, the 
action variable 



1 , 2V2^e^/^ 

Jz^ - Pzdz = . 

TT Jq 3neF 

In turn, the energy expressed in terms of Jz is 

e = [SirJzeF/VS^]'^^^ ■ 

From Eq. (jC2[) one can define the angle variable 

de ireF 
Oz = UJzt = — — t — , t 
dJz V2m*£ 



(CI) 



(C2) 



(C3) 



^t~^t', 0<t<T/2. (C5) 
TO* Zm* 

By dint of equation t — Oz/ujz, it is useful to establish 
the relation between the coordinate z and the variable 9-, 



, Q<ez<TT. (C6) 



TT \ TT 



For the triangular potential one has the semiclassical 
quantization condition Jz = hijiz + 3/4), where Uz = 
0, 1, 2, ... and the Maslov index = 3/4. Redefining 
Uz — I, one obtains the semiclassical energy levels from 
Eq.dCll 



{eFnf 



21 1/3 



2to* 



Cn.^[H^z-\)f'^ (C7) 



with Uz — 1,2, 
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